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1 Introduction 


This report investigates the use of active control to attenuate structural vibrations of 
the NASA Langley Phase Zero Evolutionary Structure due to external disturbance ex- 
citations. Hoc and structured singular value (//) based control techniques are used to 
analyze and synthesize control laws for the NASA Langley Controls-Structures Interac- 
tion (CSI) “Evolutionary Model” (hereon denoted by CEM) [Bel91,LimB]. The CEM 
structure experiment provides an excellent test bed to address control design issues for 
laxge space structures. Specifically, control design for structures with numerous lightly 
damped, coupled flexible modes, collocated and noncollocated sensors and actuators 
and stringent performance specifications. The performance objectives are to attenuate 
the vibration of the structure due to external disturbances, and minimize the actuator 
control force. 

The control design problem formulation for the CEM Stucture uses a mathematical 
model developed and refined by NASA researchers with finite element techniques. A 
reduced order state space model for the control design model is formulated from the 
finite element model. It is noted that there are significant variations between the design 
model and the experimentally derived transfer function data. These model inaccuracies 
or uncertainties take the form of unknown parameter values, natural frequencies, damp- 
ing levels, and unmodeled high frequency dynamics. Control laws are designed which 
accounted for these uncertainties and it is seen that the incorporation of errors between 
the physical system and its mathematical models into the design process is essential to 
achieve the desired performance objectives. 


2 CSI Evolutionary Model 

The structure shown in figure 1 is the CEM structure located at NASA Langley Research 
Center. The main section of the truss is approximately 50 feet long to which are attached 
two vertical appendages. Mounted on the short appendage is a 16 ft. diameter reflector 
with a circular mirror attached to its center. A laser beam whose source is attached to 
the end of a vertical appendage is pointed at the mirror. The laser beam reflection from 
the mirror is detected by a spatially fixed optical detector array (target plant) mounted 
approximately 60 ft. above the structure. 

A finite element model of the CEM structure was developed by NASA Langley re- 
searchers. The model contained 86 modes corresponding to all of the CEM structural 
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Figure 1: Schematic of the NASA Langley Controls-Structures Interaction (CSI) Evo- 
lutionary Structure (CEM) 
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Mode Frequency (rad/s) 

Mode Frequency (rad/s) 

1 

0.924 

2 

0.937 

3 

0.975 

4 

4.587 

5 

4.699 

6 

5.491 

7 

9.258 

8 

10.921 

9 

11.831 

10 

14.460 

11 

17.835 

12 

25.225 

13 

25.335 

14 

26.425 

15 

27.595 

16 

34.566 

17 

38.827 

18 

39.149 

19 

40.658 

20 

41.908 

21 

46.322 

22 

52.108 

23 

56.337 

24 

78.451 

25 

105.89 




Table T. CEM Modes and Natural Frequencies 


modes under 50 Hz. A 25 mode reduced order model of the structure, based on the 
86 mode finite element model, is obtained through a controllability and observability 
analysis. These 25 modes correspond to the set of modes below 20 Hz that significantly 
contribute to the measured accelerations and LOS displacements. The first six lowest 
frequency modes are the pseudo rigid body modes (pendulum modes) which range from 
0.14 Hz to 0.87 Hz, the seventh and eight modes are the first horizontal and vertical 
plane bending modes respectively. Table 1 contains the natural frequency of each mode 
in the CEM reduced order control design model. 

For this study, the eight air actuators are used for control and the eight nearly collocated 
accelerometers are the measurements used for feedback (see Fig. 1). A 50 th order, 8 
input and 8 output, state-space model of the CEM structure, denoted as mod25, is 
constructed from the 25 mode reduced order finite element model for control analysis. 

Model of the air actuators are taken as a constant gain in all the problem formulations 
since the bandwidth of the actuators are 50 Hz, which is significantly higher than design 
model. Similarly, the accelerometers, which have a bandwidth of 200 Hz, are treated as 
constant gains in the control designs. For further details on the evolutionary structure 
one is referred to References [Bel91]. 
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3 Control Design Models for the CEM Structure 


A measure of the accuracy of the control design model or alternatively a measure of the 
model error is needed in the control analysis and synthesis process. Transfer functions 
between the CEM control actuators and sensors were derived experimentally via Fourier 
Transform techniques by NASA researchers to access the accuracy of the control design 
model. Significant differences were noted between the experimental transfer function 
data and the 50 th state-space model, mod25, above 4 Hz. Figures Al, A2, A3, A4, A5 
and A6 in Appendix A contain plots of the experimentally derived transfer functions and 
mod25 transfer functions between air actuator 2 and acceleroemeter 2, air actuator 2 and 
accelerometer 4, air actuator 3 and accelerometer 3, air actuator 4 and accelerometer 4, 
air actuator 6 and accelerometer 6, and air actuator 8 and accelerometer 6. 

The variation between the experimentally derived transfer functions and the design 
model are accounted for in the control design process via unstructured uncertainty mod- 
els. An additive uncertainty model is used to account for neglected structural modes in 
the design model. An actuator input uncertainty model is used to incorporate any errors 
between the experimental data and the 50th state-space model which occur within the 
control bandwidth. Alternatively, these modeling inaccuracies could have been described 
by parametric uncertainty in the natural frequencies, damping values and mode shapes 
of the control design state-space model. Initially, the uncertainty models in the control 
problem formulation are restricted to unstructured multiplicative input uncertainty and 
additive uncertainty. Describing the modeling errors as parametric uncertainty in the 
natural frequencies, damping values and mode shapes of the control design model is ad- 
dressed in the paper on design of controllers to achieve line-of-sight (LOS) performance 
requirements for the CEM structure [LimB]. 

The additive uncertainty weight is selected such that the magnitude of the its transfer 
function description bounds the peaks of the neglected structural modes. Hence, the 
additive uncertainty weight is highly dependent on the level of damping assumed in the 
design model. The damping values associated with the structural modes of the CEM are 
difficult to predict, therefore based on estimates by engineers at NASA Langley, a value 
of 0.5% damping is assumed for all of the flexible modes in the control design model. 
The additive uncertainty weights are chosen accordingly. 

Eight additive uncertainty weights are formulated based on the variation between the 
control design model and the experimentally derived transfer function data. Each weight 
corresponds to the input/output relationship between an air actuator and the eight 
accelerometer measurements. The additive uncertainty weights are varied in the design 
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process to determine their effect on the closed-loop robustness and performance of the 
controllers synthesized with the //-framework. 

Three models of the CEM structure are used in the control design problem formulation. 
All of these models have 8 actuator inputs and 8 sensor outputs. mod25r is a 22 state 
model of the CEM structure, containing the first eleven structural modes between 0 and 
4 Hz. The 14 modes between 4 and 10 Hz, which are in the 50 state model mod25 but 
not mod25r, are treated as unmodeled dynamics along with all of the higher frequency 
structural modes not included in the 50 th order design model. mod25r is the design 
model used in the control problem formulations 1 through 4. 

Design model mod25r2 consists of the first eleven modes of the CEM structure along 
with the next ten structural modes, for a total of 42 states. This model describes 
the CEM structure between 0 and 6.8 Hz. As before, the neglect structural modes 
and higher frequency modes are accounted for in the control problem formulation via 
additive uncertainty models. The goal of this model is to increase the bandwidth of the 
controller to attenuate vibrational modes of the CEM structure up to 6.8 Hz. mod25r2 
is the model of the CEM structure used in control problem formulations 6 through 8. 
Design model mod25z consists of the first 15 modes of the CEM structure. The first 
eleven modes are the same as mod25r with modes 12 through 15 at 25.2, 25.3, 26.4 and 
27.6 rad/ sec. The additional modes of the structure are included near the bandwidth of 
the closed-loop system to evaluate their effect on the robustness properties of the control 
design. Control problem formulation 9 uses mod25z in the control design process. 

All of the controllers synthesized with the three models of the CEM structure are simu- 
lated with the 50 th order state-space model and implemented on the real CEM structure. 


4 The Structured Singular Value (/it) Framework 

The structured singular value (//) framework is used to analyzes the robustness of 
model sets to structured and unstructured uncertainties and synthesize controllers which 
achieve desired robust performance objectives. The model sets are constructed for the 
each control problem formulation by introducing perturbations, of the appropriate mag- 
nitude, in all uncertain parameters associated with the design model. The new design 
model encompasses a set of uncertain systems defined by a model whose parameters lie 
in the set plus the norm-bounded unstructured uncertainty. 

The //-framework is used to analyze and synthesize controller for the Evolutionary Phase 
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Zero Structure. It provides a quantitative way of measuring robustness and performance 
in the same framework. The generic model structure, figure 2, is based on linear frac- 
tional transformations (LFT). Any linear interconnection of inputs, outputs and com- 
mands along with perturbations and a controller can be viewed in this context and 
rearranged to match this diagram [Doyl]. Figure 3 (a) is a diagram of the control 
analysis problem. The analysis problem assumes a controller K has been provided and 
the bottom loop of the augmented plant P has been closed. The goal is to determine 
the stability of the closed-loop system in the presence of the norm-bounded uncertainty 
block A. The controller synthesis problem is shown in figure 3 (b). The goal is to 
design a controller K which minimizes the size of the output errors e due to a set of 
norm-bounded input v in the sense of the structured singular value, /i. 

The following are references that provide the background on /i- analysis and synthesis 
techniques [Doyl, Doy2, Pack, Mutools, Balas, BalDoyl, BalDoy2, DoyLP], 



Figure 2: General Interconnection Structure 




Figure 3: (a) Analysis and (b) Synthesis Problem 
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5 Control Objectives 


The performance objective is to attenuate the sensed vibration at accelerometer locations 
1 through 8. These sensors are nearly collocated with air actuators 1 through 8. The 
disturbances used to excited the structure enter at air thruster 1, 2, 6, and 7. The 
performance criteria is defined as minimizing the ||-||oo norm of the transfer function 
from the input disturbances to accelerometers 1 through 8. In problem formulations 
1 and 3, the performance objective is to attenuate only the structural modes which 
affect sensors 1, 2, 7, and 8. The input disturbances are modeled in the control problem 
formulation as white noise inputs to the air actuators. It was decided not to use filtered 
white noise in the problem formulation to reduce the number of states in the control 
design. 

Approximation of the input signals as white noise is based on H ^ norm requirement to 
attenuate the resonant peaks of the disturbance to error transfer function. The highest 
resonance peak will receive the most attention from the controller since it corresponds 
to the largest gain from input to output. The resonance peaks of the CEM structure 
are two order of magnitude larger than the steady state gain of the transfer function at 
high frequency. Therefore, filtering the input or performance weight isn’t required since 
the high frequency gain of the system is a factor of 100 lower than the resonant peaks. 
To verify this control design 6 through 8 filter the performance weight, which has little 
effect on the performance of the controllers. 


6 Control Problem Formulation 

The control problem formulation interconnection structure is shown in figure 4. The 
problem formulation is the same for the eight control designs with the additive uncer- 
tainty weight, Addunc, input multiplicative uncertainty weight, actu, and performance 
weight perfwt varying with each design. The goal of the control design is to minimize 
the low frequency vibration of the accelerometers due to input disturbances entering at 
air actuators 1, 2, 6 and 7. The controllers are tested by inputting a disturbance signal 
into the air thrusters for 9 seconds, terminating the disturbance signal and allowing the 
CEM structure to vibrate in the open-loop configuration for one second, on implementing 
the controller which uses the air thrusters to attenuate the induced vibrations. 

Four input disturbances are used to excite the CEM structure experimentally. Distur- 
bance signal 1, 2sin(9.24t), is input to air thruster 1 to excite the first X-Y bending 
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mode, disturbance signal 2, 0.5 sin(10.93£) is input to actuator 2 to excite the first X-Z 
bending mode, disturbance signal 3, 2sin(.974t), is input to actuator 7 to excite the yaw 
pendulum mode and the fourth disturbance signal, 0.5sin(.924£), is input to actuator 
8 to excite the X pendulum mode. The disturbance signals are treated as pure white 
noise inputs into the system. 

The uncertainty descriptions in the control design problem formulation defined the set of 
plant models in which the actual CEM structure is assumed. The performance weights 
are used to define the vibration attenuation requirements at each accelerometer location. 



Figure 4: CEM Control Problem Formulation 


7 Uncertainty Models 

Two types of uncertainty are used to account for modeling errors between the design 
model and experimentally derived transfer functions. An input multiplicative uncer- 
tainty accounts for variations between the model and the real structure within the con- 
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trol bandwidth and an additive uncertainty model accounts for unmodeled dynamics. 
The additive uncertainty ensures that the controller will not destablize the unmodeled 
modes of the structure, provided the nominal structural model and the additive uncer- 
tainty description accurately represent the dynamics of the CEM structure. 

The input multiplicative uncertainty weight is selected to be a constant in first seven 
of the eight control designs. The reasoning behind this is that the input uncertainty 
weight is only important within the control bandwidth. Outside the control bandwidth 
the additive uncertainty weight dominates the control problem and it is difficult to 
determine the differences in the modeling error of the pendulum modes as compared 
with the first bending modes, hence the same level of uncertainty is selected for both. 
Another approach to account for the modeling error is to include parametric uncertainty 
in the natural frequency and damping coefficients in the control design model. This is the 
approach taken to design controllers to minimize LOS errors due to external disturbances 
excitations [LimB], 

Experimental results have indicated that as long as some input uncertainty is included 
in the problem formulation, the controllers will be robust to errors within and the perfor- 
mance will not be significantly degraded [Balas, BalDoy2]. Determining the correlation 
between a minimal level of uncertainty and a large amount of uncertainty is an area of 
research which will be pursued this year by the principal investigator. 

Table 2 contains the actuators, sensors, uncertainty weights and performance weights 
used in the control problem formulations to synthesize the eight controllers. Note that 
there is no control design 5. \i - synthesis techniques based on D — K iteration are used to 
design the controllers to achieve robust performance for the given problem formulation. 
The eight controller designed all had on the order of 80 states prior to being reduced. 
The large number of states is due to the order of the conrol design model, additive 
uncertainty weights and D-scaling used in the D — K iteration process. The reduced 
order controllers for each design, the controller order reduction is done using balanced 
realization techniques, are listed in Table 2. 


8 Control Design Results 

All eight of the controllers designed using the ^/-synthesis techniques resulted in a stable 
closed-loop system and performed very well when implemented on the experimental CEM 
structure. Control design 9 achieved the best performance of any controller previously 
tested. These results indicate that the /i-framework is ideally suited for the analysis 
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Design 

Model 

Sensors 

Actuators 

Input 

Uncertainty 

Additive 

Uncertainty 

Performance 

Weight 

Controller 

Order 

1 

mod25r 

1,2, 7, 8 

1,2, 7, 8 

10% 

addwt 1,2, 7, 8 

10:1 

28 

2 

mod25r 

1 -> 8 

1 — ► 8 

10% 

addwt 1 — ► 8 

10:1 

36 

3 

mod25r 

1,2, 7, 8 

1,2, 7, 8 

2% 

addwt 1,2, 7, 8 

10:1 

30 

4 

mod25r 

1 8 

1 -► 8 

1% 

addwt 1 — > 8 

10:1 

40 

6 

mod25r2 

1 

00 

1 8 

1% 

addwt 1 — » 8 

10:1, filters 

50 


mod25r2 

1 -> 8 

1 -» 8 

10% 

addwt 1 — > 8 

8:1, filters 

50 


mod25r2 

1 -»■ 8 

1 8 

100 (s+20) 

s+2000 

addwt 1 — > 8 

8:1, filters 

38 


mod25z 

00 

T 

1 -4 8 

10% 

addwt 1 — > 8 

perfwt9 

44 


Table 2: Control Problem Formulation Uncertainty and Performance Weights 


and synthesis of controllers for flexible structures with numerous, lightly damped modes 
and multiple actuators and sensors. Experimental results of open-loop system and the 
eight controllers are located in Appendix A, figures A7 to A42. The experiments on the 
CEM structure use four disturbance signals input to actuator 1, 2, 7 and 8 to excite the 
structure. The open-loop accelerometer responses of sensor 1 through 8 and actuator 1 
through 8 inputs are shown in figures A7 through A10. 

Controller 1 was designed using only four of the eight air actuators available. It was 
implemented on the actual CEM structure and performed as well as the previous best 
controller. Controller 2, used the same uncertainty and performance weights as control 
design 1, employed all eight actuators and sensors for control. The improvement in 
performance was limited, although only one set of disturbance excitations were used in 
the experiments. The most notable improvement with controller 2 was the attenuation 
of vibration at accelerometers 3 and 6. This is to be expected since accelerometers 3 and 
6 are almost collocated with air actuators 3 and 6, which were used in control design 2 
but not in control design 1. 

The role the input multiplicative uncertainty level plays is evident when comparing 
the performance of controller 1 with that of controller 3. Controller 3 is designed the 
same as controller 1 except that the input multiplicative weight is selected to be 2% 
uncertainty compared with 10% uncertainty. The performance of controller 3 at all of 
the accelerometer locations is slightly improved. The increase in performance comes at 
the expense of increased actuator force, most notably from air thrustor 1 and 8. Results 
from the implementation of controllers 1 and 3 indicate that the uncertainty models 
used in the problem formulations accurately represent the physical CEM structure. The 
theoretical results for each design parallel the experimental data exactly. 
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Controller 4 uses the same problem formulation as control design 2 with the input 
multiplicative uncertainty weight choosen to be 1% as opposed to 10%. Similar results 
are obtained to controller 1 and 3. Controller 4 has increased the attenuate of low 
frequency virbation at the accelerometer locations at the expense of increased actuator 
control forces. It is interesting to note the high frequency content in the accelerometer 
3 signal. The increased controller bandwidth in design 4 causes several high frequency 
modes to be excited. 

The problem formulation used to synthesize controller 6 uses a higher order design model, 
different additive uncertainty models reduced in magnitude from designs 1 through 4 and 
first order filters on the performance weights. The performance of controller 6 is similar 
to controller 1, though controller 6 has significantly more high frequency content in the 
accelerometer signals. The additive uncertainty models used in design 6, allowed the 
controller to have a higher bandwidth leading to increased high frequency accelerome- 
ter signals. The higher bandwidth controller lead to poor robustness characteristics for 
controller 6. Similar problems were noted in control design 7. The input multiplicative 
uncertainty problem formulation was modified to 10% to account for the change in the 
additive uncertainty models. The results indicate that accounting for unmodeled dy- 
namics via input multiplicative is not a good approach to the control design problem 
for lightly damped flexible structures. Based on intuition, one would expect that a very 
large multiplicative uncertainty weight would be required on the actuator in the neigh- 
borhood of an unmodeded flexible mode of the structure. A level of 10% multiplicative 
uncertainty would be inadequate to assure stability of the unmodeled flexible modes. 
Controller 7 appears to have destabilized a structural mode of the system around 7 Hz. 

Control design 9 represents the controller with the highest level of performance. A high 
order model of the structure is used in the design process, with the goal to attenuate 
the structural mode at 7 Hz. The performance weight perfwt is modified to reflect this 
goal and the additive uncertainty weight is modified accordingly. The controller atten- 
uates the low frequency modes of the CEM structure quickly with little high frequency 
content to the acceleroemeter signals. Unfortunately, due to the limited accuracy of the 
control design model above 4 Hz, the modeling errors between the design model and 
experimentally derived transfer function data can be seen in Appendix A figures Al 
through A6, controller 9 is unable to attenuate the structural mode at 7 Hz. In spite of 
this shortcoming, controller 9 achieved the high level of performance of any controller 
implemented on the experimental CEM structure. 

It is of interest to noted that all of the controllers implemented were obtained using 
balanced realization model reduction techiques. Although the original controllers were 
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on the order of 80 states, the controllers implemented were reduce to approximately 40 
states prior to implementation. This indicates that a number of states are not necessary 
in the control designs. It would be of interest to invesitgate how the controller order 
effect the achievable performance given the same problem formulation. 


9 Summary 

Hoo and structured singular value (/ 1 ) based control design methods are used to ana- 
lyze and synthesize control laws for the NASA Langley Controls-Structures Interaction 
(CSI) “Evolutionary Model.” The CEM structure had numerous lightly damped, cou- 
pled flexible modes, collocated and noncollocated sensors and actuators and stringent 
performance specifications which lead to a difficult control problem. The control design 
problem was further complicated by the model errors noted between the model used for 
control design and experimentally derived transfer function data. These model inaccu- 
racies or uncertainties are accounted for in the control design process via unstructured 
uncertainty models in addition to the performance specifications on the attenuation 
of structural vibration due to external disturbances. Eight controllers are synthesized 
for a variety of input multiplicative uncertainty, additive uncertainty and performance 
weights. 

The controllers exhibited excellent correlation between the theoretically predicted and 
experimentally achieved results especially associated with the low frequency modes of 
the CEM structure. This was expected due to the error between the design model and 
experimenta at frequenies above AHz. The /i- synthesis techniques resulted in controllers 
which traded off between robustness and performance given the specific problem formu- 
lation and achieved outstanding performance when implemented on the CEM structure. 
It was seen that unstructured uncertainty models of the modeling error were more than 
adequate in the control problem formulation. 


10 Research in Progress 

Parameteric uncertainty models of variation in natural frequency and damping levels is 
explored in detail by Dr. Kyong Lim at NASA Langley and myself in the paper [LimB]. 
The same and structured singular value techniques are applied to the CEM structure 
to minimize the LOS error due to external disturbances. It was found that because the 
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CEM structure is very lightly damped, the robust performance for fine pointing was 
highly dependent on the accuracy of the structural frequencies. This was not unexpected 
since small variations in natural frequencies lead to large gain and phase variations in the 
neighborhood of those natural frequncies in lightly damped systems. These results were 
not verified on the experimental structure due to the controllers being synthesized after 
the CEM structure was disassembled. Modeling of system with parametric uncertainty 
in structural natural frequencies and damping values is currently being pursued by Dr. 
Kyong Lim and myself and will be an main area of research in the coming year. 

Two Ph.D. graduate students at the University of Minnesota, Richard Lind and Arun 
Kumar are being support by this NASA grant. They are pursuing research topics appli- 
cable to the area of control of flexible space structures. Mr. Lind has developed realtime 
control software for a DSP processor which is housed within a Macintosh FX. The Mac- 
intosh has 16 A/D channels and 8 D/A channels for realtime control. The realtime 
control can be implemented at 2 kHz while input and output channels are graphically 
on the Macintosh screen. A two mass, three spring experiment has been developed to 
test out the realtime software. This experiment is being used by Mr. Kumar to inves- 
tigate the tradeoff between collocated and noncollated control. Mr. Lind is also doing 
theoretical research on calculating optimal constant D scalings for the p, upper bound 
via convex optimization. This work may have a significant impact in the synthesis of 
optimal (i controllers for the full information control design problem. 

Mr. Kumar is investigating the tradeoff between collocated and noncollocated control 
design. The idea is to develop a quantitative relationship between the degree of un- 
certainty in a system and the need for collocated control for robustness purposes. Mr. 
Kumar has spent a majority of this concentrating on coursework will be spending the 
spring quarter primarily doing research. 

We are in the process of building a large flexible structure experiment in the Aerospace 
Controls Lab at the University of Minnesota. The structure will investigate modeling, 
system identification, model validation and control analysis and synthesis issues for 
flexible structures. The structure will have a slewing component which will introduce 
nonlinearities into the problem formulation. This experiment will be used to investigate 
control design issues for flexible structures prior to testing these ideas on the CEM 
structure at NASA Langley. 
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Experimental Response of Accel 5: Open-loo 
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Experimental Response of Accel 5: Controller kdesl 
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Figure A12 
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Experimental Response of Accel 2: Controller kdes2 
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Experimental Response of Accel 3: Controller kdes2 
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Experimental Response of Accel 4: Controller kdes2 
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Figure A15 
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Experimental Response of Accel 5: Controller kdes2 



Time (seconds) 


Experimental Response of Accel 7: Controller kdes2 
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Experimental Response of Accel 8: Controller kdes2 
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Experimental Response of Accel 5: Controller kdes3 
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Figure A20 
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Experimental Response of Accel 7: Controller kdes4 
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Figure A24 
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Figure A25 
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Experimental Response of Accel 5: Controller kdes6r 
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Experimental Response of Accel 5: Controller kdes9r 
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Figure A36 
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